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Abstract 

In this study, interacting crack growth in an infinite plate is 
analyzed with new, fast and accurate Boundary Cracklet 
Method (BCM) developed by Phoenix and Yavuz. An 
interior crack is under consideration to examine its 
propagation because of cyclic loading which is very common 
for aerospace, naval and civil engineering structures. BCM is 
used to determine the overall stress field as well as stress 
intensity factors for crack tips and singular wedges at crack 
kinks. BCM uses integral equations expressed in terms of 
unknown edge dislocation distributions along crack lines. 
These distributions derive from an accurate representation 
of the crack opening displacements using power series basis 
terms obtained from wedge eigenvalue analysis, which leads 
to both polynomial and non-polynomial power series. The 
process is to choose terms of the series and their exponents 
such that the tractions on the crack faces are virtually zero 
compared to the far field loading. Applying the method 
leads to a set of linear algebraic equations to solve for the 
unknown weighting coefficients for the power series basis 
terms to make no use of numerical integrations. A simple 
crack growth emanating from a triangular hole in an infinite 
plate is analyzed. The fatigue crack growth is assumed to 
follow Paris-Erdogan Law. The results are compared to 
those of other numerical methods. A parametric study is 
performed via graphs and tables to demonstrate the ability 
of BCM in analysis of fatigue crack growth. 
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Introduction 

Interaction and propagation of cracks is a significant 



concept in the design and analysis of aerospace, naval, 
and civil structures. Sharp corners at the holes bring 
about stress localizations that reduce strength and 
fatigue life. 

In the current study, a simple straight crack growth 
emanating from a triangular hole in an infinite plate is 
analyzed with new, fast and accurate Boundary 
Cracklet Method (BCM) developed by Phoenix and 
Yavuz (Yavuz et al. (2006-2012)). An edge crack 
around the triangular hole is under consideration to 
observe its propagation. The fatigue crack growth is 
assumed to follow Paris-Erdogan Law. Triangular hole 
is modeled by straight cracks. Thus, they form an 
array of cracks within the structure. The study is 
performed via graphs and tables to demonstrate the 
ability of BCM in analysis of fatigue crack growth. 

Introductory Remarks on Boundary Cracklet 
Method (BCM) 

BCM is very useful to determine the overall stress field 
as well as stress intensity factors for crack tips and 
singular wedges at crack kinks. BCM uses integral 
equations expressed in terms of unknown edge 
dislocation distributions along crack lines. These 
distributions derive from an accurate representation of 
the crack opening displacements using power series 
basis terms obtained from wedge eigenvalue analysis, 
which leads to both polynomial and non-polynomial 
power series. The process is to choose terms of the 
series and their exponents such that the tractions on 
the crack faces are virtually zero compared to the far 
field loading. Applying the method leads to a set of 
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linear algebraic equations to solve for the unknown 
weighting coefficients for the power series basis terms 
to make no use of numerical integrations unlike in 
other methods. That's why solution takes just a few 
seconds on a PC. 

In the application, all boundaries and internal cracks 
are represented by cracklets. There are three types of 
cracklets (CI, C2 and C3) which have their own typical 
opening profile approximation functions. In the 
formulation, b is the opening profile approximation 
function, P's are polynomial subf unction family, W's 
are wedge subfunction family formed by rational 
powers of the position (t) in each cracklet. These 
rational powers are eigenvalues (q) calculated from 
Williams Wedge Analysis. Depending on the wedge 
angle, there may be only one (in normal mode, mode I) 
or two (in both normal and tangential directions, 
mode I and mode II) eigenvalues. These eigenvalues 
are approximated by a ratio to evaluate all integrals 
analytically which makes the algorithm fast. These 
subfunctions are given for two modes of crack 
openings: i=l for mode I, i=2 for mode II. All details of 
these opening profile functions and analytical and 
computational aspects of the method can be found in 
Yavuz et al (2006). 

These three typical cracklets and corresponding 
approximation functions of opening displacement 
profiles are shown in Fig. 1. Cracklet type CI is a 
cracklet with tip ends: A single straight crack is 
represented by this type of cracklet. Here only two 
wedge subfunctions with the power of (Vi) are used 
for each tip end points. 

bf (t) = W,f (t la t ) + (1-t/a,) (1) 

Cracklet type C2 is a cracklet with one wedge end and 
one tip end: A V-crack has two this type of cracklets. 
At tip end point B, the wedge eigenvalue is (Vi) so we 
have only one wedge subfunction family and for the 
other wedge end point A, there are one polynomial 
and one or two wedge subfunction families as C2. 



(0 = 



PT\t 1(1,) + WT (t laM +WT(t I a. ) + W° 

] Pi J 1 \ Pi J 1 1 



\l-tla t ) (2) 



Cracklet type C3 is a cracklet with only wedge ends: A 
kinked crack is represented by three cracklets and the 
yellow one is this type of cracklet. The opening profile 
approximation (b) may have from three to five 
subfunction families depending on the wedge angle at 
end points A and B. 

b f it ) = P f B W it /a l )+W^f(t/a l )(+ W^fit/aOyW^fil-t/a^l+W^Jil-t/a^ 

(3) 



Even problems with curved boundaries and cracks can 
be solved by using small cracklets representing curved 
shapes. Since a small C2 type cracklet is added to each 
present crack's end to capture fatigue crack 
propagation after specific number of cyclic loads, all 
calculated SIFs at the end points become SIFs of the 
whole propagated crack configuration. 




(a) CI 




(b) C2 




n h yi 

(c) C3 

FIG. 1 CRACLET TYPES 

Crack Growth and Propagation Model 

Fatigue life calculations and crack growth behaviour 
are determined by Paris-Erdogan Law (Paris and 
Erdogan (1963)). Paris-Erdogan Law can be expressed 
by the formula given below, 



dc 
dN 



■■C(AK)" 



AK = K„ 



(4) 



(5) 



where dc/dN represents crack growth rate, AK is stress 
intensity range, C and m are material constants. 
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Constant amplitude cyclic loading with Omin=0 is 
considered. Loading condition is illustrated in Fig. 2. 
Here, R is the stress ratio, o m is the midrange stress, 
and a a is the stress amplitude. 



R- 



(6) 
(7) 



Constant Amplitude Loading 



Stress 




Time 



FIG. 2 CONSTANT AMPLITUDE LOADING WITH R=0 AND 

Omin=0 

Crack propagation and fatigue life estimation analysis 
is performed for the material parameters given in 
Table 1 (Yan (2006)). 

TABLE 1 MATERIAL PARAMETERS 



Material 
Parameters 


Definition 


Values 


G 


Shear Modulus [MPa] 


26918.64 


N 


Poisson's Ratio 


0.321 


Kic 


Fracture toughness [MPa.mml/2] 


1137.96 


AKth 


SIF threshold [MPa.mml/2] 





C 


Material constant [mm/cycle] 


1.039xl0 10 


m 


Material constant 


2.7438 



Two different stress ranges (Act) are considered for the 
fatigue life and crack propagation analysis. Infinite 
plate is subjected to uniaxial tension in plane stress 
conditions. Initial crack orientations are taken as 0°, 30°, 
45°, and 60°, respectively. For a 0° oriented crack 
subjected to uniaxial tension, Mode I condition is 
dominant. Since the orientations of the cracks that 
form triangular hole are different from 0°, Mode II 
effects must be included. It is assumed that the fatigue 
crack growth is governed by the formula given below 
(Sih and Barthelemy (1980)), 



dc 
dN 



:B(AK eff ) S 



(8) 



where 



AK eff ^-cos^-^.a + cos^^-SA^sin^] (9a) 



B=C, s=m (9b) 

where Kejf is the effective stress intensity factor, 9o is 
the crack growth angle, B and s are material constants. 
Crack propagation is performed by setting a constant 
crack propagation length for each step. Crack 
propagation angle is determined by considering the 
direction of the maximum tangential stress around the 
crack tip (Liang et al. (2003)). 

c. +1 = c t + Ac f (Ac f : constant) {i = 1,2,3,..., J) (10) 

f \ 
-2 



(0 o ) ; =2tan- 




(11) 



Fatigue life calculation is Followed by the formula 
given below, 



N 4 = 



B(AK eff \ 



(* = 1,2,3,...,/) (12) 



Case Study (An Infinite Plate with a 
Triangular Hole Subjected to Uniaxial 
Tension) 

The loading conditions and the geometry of the 
problem are shown in Fig. 3. 

The stress intensity factors are compared with the 
results found in the literature to show the accuracy of 
BCM (Table 2). The results are found to be in a good 
agreement with existing solutions in studies 
Murakami (1987) and Yan (2005). The crack 
propagation is analyzed for four different initial crack 
orientations and two different stress ranges. The 
parameters are given in Table 3. The shape of the crack 
emanating from a triangular hole and its propagation 
is shown in Fig. 4. As expected, the crack becomes 
perpendicular to the applied load direction regardless 
of its initial orientation. The effective stress intensity 
factor versus crack length is given in Fig. 5 and Fig. 6. 
The higher stress amplitudes cause the higher stress 
intensity factors. The stress intensity factors are found 
to be higher for the lower values of initial crack 
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orientations. In the case of lower values of initial crack 
orientations, it is found that the number of cycles to 
reach a given crack size is lower compared to the one 
for the higher values of initial crack orientations (Fig. 
7). 



Crack Propagation - Triangular Hole 





FIG. 4 CRACK PROPAGATION GRAPH 



\* \? \* ^ j^ 1 



FIG. 3 AN INFINITE PLATE WITH A TRIANGULAR HOLE 
SUBJECTED TO UNIAXIAL TENSION 

TABLE 2 COMPARISON OF THE STRESS INTENSITY FACTORS 
FOR A TRIANGULAR HOLE WITH CRACK (6=0°) 



Crack Length vs. AK gff Ao=1 50.42 MPa 





Murakami (1987) 


Yan (2005) 


BCM (current) 


cla 


Ki/(nc) - 5 


Ki/(nc) - 5 


Ki/(nc) - 5 


0.05 




3.90 


3.9148 


0.1 


2.76 


2.8348 


2.8414 


0.2 


2.064 


2.0650 


2.0686 


0.4 


1.525 


1.5238 


1.5266 


0.6 


1.297 


1.2959 


1.2975 


0.8 


1.168 


1.1671 


1.1688 


1.0 


1.085 


1.0786 


1.0863 



TABLE 3 PARAMETERS OF THE CASE 





Case 


Definition 


Ao 


150.42, 300.83 


Stress ranges [MPa] 





0°, 30°, 45°, 60° 


Orientation of the initial 
crack 


a/a 


0.2 


Ratio of the length of the 
initial crack to the edge 
length (i=0) 


Ada 


0.1 


Ratio of the crack length 
increment to the initial crack 
length (i=0) 



'E 



300 
250 
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2 2.5 
Crack Length (mm) 



FIG. 5 CRACK LENGTH VS AKeff (zla=150.42 MPa) 

Crack Length vs. AK gff Ao=300.83 MPa 
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2 2.5 
Crack Length (mm) 



FIG. 6 CRACK LENGTH VS AKeff (zlc7=300.83 MPa) 
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Crack Propagation vs. Fatigue Life 
3.5 I , 1 , , 




Number of Cycles 

FIG. 7 CRACK LENGTH VS NUMBER OF CYCLES 

Conclusion 

The Boundary Cracklet Method is used to anticipate 
the stress intensity factors for the cracks emanating 
from a triangular hole during the cyclic loading. The 
stress intensity factors calculated for several initial 
crack orientations are found to be in an agreement 
with the ones found in the literature. The fatigue crack 
growth is assumed to follow Paris-Erdogan Law. The 
results are obtained for the propagation of cracks 
emanating from a triangular hole. As the number of 
cycles increases, the crack becomes perpendicular to 
the applied load direction regardless of its initial 
orientation as expected. Also, the higher stress 
amplitudes cause the higher stress intensity factors. 
The stress intensity factors are found to be higher for 
the lower values of initial crack orientations. In the 
case of lower values of initial crack orientations, it is 
found that the number of cycles to reach a given crack 
size is lower compared to the one for the higher values 
of initial crack orientations. The method presented 
here predicates the crack propagation for the cracks 
emanating from a triangular hole very well. 
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